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ABSTRACT 

Most of gamma-ray bursts (GRBs) models have predicted that the intrinsic 
isotropic energy is limited to below ~ 10 53 ~ 54 ergs. Recently claimed high red- 
shifts, correlation with supernovae, and connection to the cosmic star formation 
activity point to a different energy requirement and source evolution history possi- 
bly with strong beaming. We study the effects of the beaming-induced luminosity 
function on statistics of observed GRBs, assuming the cosmological scenario. We 
select and divide the BATSE 4B data into 588 long bursts (T 90 > 2.5 sec) and 
149 short bursts (T 90 < 2.5 sec), and compare the statistics calculated in each 
subgroup. The (V/V max ) of the long bursts is 0.2901 ± 0.0113, and that of the 
short bursts is 0.4178 ± 0.0239. For luminosity function models, we consider a 
cylindrical-beam and a conic-beam. We take into account the spatial distribution 
of GRB sources as well. A broad luminosity function is naturally produced when 
one introduces beaming of GRBs. We calculate the maximum detectable redshift 
of GRBs, £ max - The estimated z max for the cylindrical-beam case is as high as 
~ 14 (a = 1.0) and ~ 6 (a = 2.0) for the long bursts and ~ 3 (a = 1.0) and ~ 1.6 
(a = 2.0) for the short bursts where a is the photon index. The large z max value 
for the short bursts is rather surprising in that the (V/V max ) for this subgroup 
is close to the so-called Euclidean value, 0.5. We calculate the fraction of bursts 
whose redshifts are larger than a certain redshift z', i.e. f >z '. When we take 
z' = 3.42 and apply the luminosity function derived for the cylindrical-beam, the 
expected f >z > is ~ 75 % (a = 1.0) and ~ 50 % (a = 2.0) for long bursts. When 
we increase the opening angle of the conic beam to A6 = 3°.0, f >z > decreases 
to ~ 20 % (a = 1.0) at z' = 3.42. If we assume a = 2.0, the conic-beam with 
A9 = 3°.0 can not explain the redshift distribution of the observed GRBs. We 
conclude that the beaming- induced luminosity functions are compatible with the 
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redshift distribution of the observed GRBs although the apparent "Euclidean" 
value of (V/V max ) might be explained by the standard model. 

Subject headings: cosmology: theory - gamma rays: bursts - methods: statistical 

1. Introduction 

The observation of high-redshift GRBs such as GRB970508, GRB971214, and 
GRB980703 (Metzger et al. 1997; Kulkarni et al. 1998; Djorgovski et al. 1998) confirmed 
the cosmological origin of GRBs (Paczyhski 1986; Meegan et al. 1992). High redshifts 
(z > 1) of GRBs would imply a much larger apparent energy budget than the conventional 
energy budget ~ xo 53-54 ergs unless some significant beaming is involved (Kumar 1999; 
Kumar and Piran 2000). The isotropic energies of observed GRBs with known redshifts 
range from ~ 10 50 ergs to ~ 10 54 ergs. A broad luminosity function can be a reasonable 
answer for this wide range, and can be naturally expected if we assume beaming of GRBs. 
The suggested correlation between GRBs and certain types of supernovae (Ib/c) implies 
that the total energy available for GRBs could be rather limited at least in the low energy 
GRBs L int ~ 10 49 erg/sec (Wang and Wheeler 1998; Woosley et al. 1999, see also Kippen 
et al. 1998; Graziani et al. 1999). If we consider the supernova explosion as a universal 
GRB source, then the reported GRB-supernova correlation implies that the high redshift 
GRBs must be strongly beamed with the beam opening solid angle, which is approximately 
given by 

dQ ~ (L int /L obs ) ~ 10- 4 (L int /10 49 erg/sec)(L obs /10 53 erg/sec)- 1 . (1) 

When we consider beaming, a broad luminosity range and a stringent constraint on 
the intrinsic energy are anticipated. Yi (1994) considered the effect of the luminosity 
function due to the cylindrical-beam with the inhomogeneous distribution of GRB sources 
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on statistics of GRBs. He analytically showed that the luminosity function follows a power 
law with the index —4/3 in the very narrow beam model when a = 1.0. His result has 
been confirmed by Mao and Yi (1994) and Chang and Yi (2000). Mao and Yi (1994) 
studied the statistical properties of luminosity functions derived using a conic beam with 
a uniform spatial distribution of GRB sources. They examined the cumulative probability 
distribution of the peak count rates. They found that the maximum redshift increases as a 
product of the Lorentz factor and the opening angle decreases. We attempt to generalize 
the luminosity function induced by the cylindrical-beam by introducing the conical shape 
in such a beam. If the progenitors of GRBs were produced in the late stage of massive stars 
(Woosley 1993; Paczyhski 1998; MacFadyen and Woosley 1999), their spatial distribution 
may follow the star-formation rate (SFR) of massive stars (Totani 1997; Wijers et al. 
1998; Blain and Natarajan 2000). We use the recently observed SFR data to derive the 
possible number density distribution of GRB sources (Steidel et al. 1999). 

In this paper, we study how beaming-induced luminosity functions affect the statistics 
of the observed GRBs in the BATSE 4B catalog (Paciesas et al. 1999) considering the 
two different spatial distributions of GRB sources: an SFR-motivated distribution and the 
uniform distribution. We employ two different SFR-motivated distributions. We discuss 
the effect of the photon index a (— 1.0 and 2.0) in the conlusion. We assume a flat universe 
with no cosmological constant, and adopt the Hubble constant H = 50km/sec/Mpc (cf. Yi 
1994; Mao and Yi 1994). We analyze 775 observed bursts in the BATSE 4B catalog in order 
to improve the statistical significance. Assuming the bimodal distribution of durations of 
GRBs (Kouveliotou et al. 1993; Lamb et al. 1993; Mao, Narayan and Piran 1994; Katz 
and Canel 1996) we divide the total sample into two subgroups according to their duration 
and study the statistics of each subgroup. 

The selection criteria applied to observational data are described in section 2. 
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Luminosity functions due to beaming are discussed in section 3. We present results and 
discussions in section 4. We summarise our conclusion in section 5. 

2. Classification of Observational Data 

We adopt the bursts in the BATSE 4B catalog (Paciesas et al. 1999) and calculate 
their (V7V max ) for selected GRBs (Schmidt et al. 1988), using fluxes in channels 2 and 3. 
We choose bursts detected in 1024 ms trigger time scale with the peak count rates satisfying 
Cmax/Cmin > 1- Applying these criteria, we select 775 bursts among 915 bursts in the 
BATSE 4B Cmax/Cmin table. The (V/V max ) value for this total sample is 0.3177 ± 0.0102. 
Then we divide the selected sample into two subgroups according to the burst durations 
(T 90 ), which is motivated by the reported bimodal structure in GRB durations (Kouveliotou 
et al. 1993; Lamb et al. 1993; Mao, Narayan and Piran 1994; Katz and Canel 1996). We 
define bursts with T 90 > 2.5 sec as long bursts, and those with T 90 < 2.5 sec as short bursts. 
Their estimated (V/V max ) values are 0.2901 ± 0.0113 (588 bursts) and 0.4178 ± 0.0239 
(149 bursts), for long bursts and short bursts respectively. The number in parentheses 
indicates the number of bursts belonging to each subgroup. The (V/V max ) values obtained 
from the two subgroups show deviations from the values which are expected for a uniform 
distribution in a Euclidean space. In our subgroups, the long bursts have a larger deviation 
from the Euclidean value. It agrees with the claim by Tavani (1998). According to Tavani 
(1998), bursts with long durations (T 90 > 2.5 sec) and hard spectra (HR 32 > 3.0) show the 
largest deviation from the so-called Euclidean value. 
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3. Beaming and Luminosity Functions 

We assume two kinds of number density distributions of GRB sources: (i) the 
SFR-motivated distribution n SFR (z) and (ii) the uniform distribution n u (z) = n o u . We 
assume a two-sided Gaussian function as usfr(z) (Kim, Yang, and Yi 1999), and fit n SFR (z) 
to the extinction-corrected SFR data (Steidel et al. 1999). The observed SFR gradually 
decreases or even keeps flat over a wide redshift range beyond a certain redshift (cf. Madau 
et al. 1998). The SFR-motivated number density distribution is approximated by 

n S Fn(z) = n 0iSF R exp - {z ~ z f , (z < z c ) (2) 



™ ,sfr exp 



Azr 

(^c) 2 

Az 2 2 



, (z > Z c ), 



where Azi and Az 2 are widths of the fitting functions for the left part and the right part 
and z c = 1.5 is a critical redshift. All parameters are determined by fitting the observational 
data. We set Azi — 1 [z < z c ) and we also effectively fix nsFR( z ) to n o S FR beyond z c . 
The proportional constants n 0)SFR and normalized to satisfy the number from the 

observational data. 



3.1. Cylindrical Beaming 

We consider the beaming-induced luminosity function following Yi (1993, 1994), 
where the relativistic beam has a constant cross-sectional area, i.e. a perfectly collimated 
cylindrical beam. Under this circumstance, a bulk motion of all the particles emitting the 
photons is parallel to the cone axis. In this case, the beam opening angle is irrelevant for the 
luminosity function. We assume that the beamed emission from "standard bursts" which 
have the same intrinsic luminosity L int . The advantage of this model is that it physically 
defines a range of the luminosity, i.e. [L min , L max ], for a given photon index a, and the 
Lorentz factor 7. An apparent luminosity function due to the uniformly distributed 
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beam in space is given by 

m = ± L % L-W (3) 

where (3 = (1 - 7" 2 ) 1/2 and p = a + 2. We adopt a = 1.0 (Mallozzi et al. 1996) and 2.0 
(Yi 1994) and the Lorentz factor 7 = 100 (see Piran 1999). If we assume a bi-polar jet, the 
minimum and the maximum luminosities are naturally given as follows: 

Lmin = £int-B m in> -^max = -^int-B max ' (4) 

where B min = 7 -1 and B max = [7(1 — (3)}~ 1 (Blandford and Konigle 1979). The intrinsic 
luminosity L int for each subgroup is determined by the (V/V max ) test. The obtained L int 
values are ~ 10 43 erg/sec for the specific set of adopted beaming parameters. This is a much 
smaller value than that required in the standard candle model, which is L ~ 10 51 erg/sec. 

We calculate two statistical quantities using the beaming-induced luminosity function 
explained in this subsection. First, we calculate the (V/V ms , x ) as a function of threshold 
flux F th (Che et al. 1999). The threshold flux is obtained when we put the maximum 
detectable redshift z max in 



F = + (5) 

\iQ1xc 2 ) [(i + ^) i /2_i]2' y°> 

where L is the luminosity, c is the speed of light, H is the Hubble constant, and a is the 
photon index. After we obtain a set of F th , we calculate the (V/V max ) for each threshold 
flux. The (V/V max ) curve converges to 0.5 as F th increases sufficiently large. We convert 
flux to the number of bursts and plot (V/V max ) as a function of the number of bursts for all 
bursts and long and short bursts in Fig. 1 (a) for a = 1.0 (thin lines) and a = 2.0 (thick 
lines). The (V/V max ) vs. number of bursts relation essentially gives the same information 
that one obtains from the conventional logN — logP plot. Our (V/V mSuX ) test does not 
include the threshold effect due to the detection efficiency. 
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Second, we calculate the fraction of GRBs. We define the fraction of bursts located at 

a redshift larger than z' as 

Number of bursts with a redshift larger than z' 
Total number of bursts 

In Fig. 1 (b), we plot the fraction of bursts as a function of redshift for each burst subsample 

for the cylindrical-beam case with different number density distributions. 



3.2. Conic Beaming 

As a more realistic example, we adopt conically-beamed emission from GRBs, which 
was studied by Mao and Yi (1994), and extended by Chang and Yi (2000). We consider 
the luminosity function produced using three different opening angles, i.e. A9 = 0°.l, 1°.0, 
and 3°.0. We calculate (V/V maiX ) and the f >z > for each subgroup with the derived luminosity 
function. The bulk Lorentz factor is set to be 100. According to Mao and Yi (1994), 
the probability that we observe the bright bursts rapidly increases as the opening angle 
increases, while relatively dim bursts are not as detectable as brighter ones. Therefore, as 
the opening angle increases, the derived luminosity function becomes similar to that of the 
standard candle case, in which all bursts have the same maximum luminosity given by this 
luminosity function. When the opening angle is small (A6 <C I/7), the luminosity function 
gives the same result as the cylindrical beaming case. Using this luminosity function, we 
calculate the same statistics as we did in previous subsection. Results are shown in Figs. 
2-5 for a = 1.0 (thin lines) and a = 2.0 (thick lines) in the cases of each sample. 



4. Results and Discussions 



The left panels in Figs. 1-5 (a) show the (V/V max ) as a function of the number of 
bursts. We show results obtained from different luminosity distributions and number 
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density distributions. The thin lines in Fig. 1-4 represent the results obtained with a — 1.0 
and thick lines a = 2.0. In Figs. 1-4, we assume that the SFR-motivated number density 
distribution, which is flat beyond a redshift larger than the critical redshift, z c . Finally, 
we show our results obtained with the SFR-motivated number density distribution, which 
gradually decreases and a = 1.0 in Fig. 5 (a). In all cases, we assume that n(z) rapidly 
increases if z < z c . The dotted lines represent the observed (V/V max ) curves with the ±3c 
bound. The solid lines indicate the theoretical results derived from our luminosity functions 
with the SFR-motivated number density distribution of GRB sources. We plot the results 
obtained for the luminosity functions with the uniform distribution of burst sources as 
dashed lines, for comparaison. Apparently all the luminosity functions we have studied 
satisfy the observed (V/V max ) curve at the 3a significance level. Although the number 
density distribution of GRB sources is believed to correlate with the SFR, we cannot rule 
out the uniform distribution of burst sources in this analysis (see also Krumholz et al. 
1998). 

The right panels in Figs. 1-5 (b) show f >z > in terms of z' . The y-axis corresponds 
to the fraction of bursts which have redshifts larger than a certain redshift z' . Solid lines 
result from the SFR-motivated number density distribution for the spatial distribution of 
GRBs and dashed lines the uniform distribution. Filled squares in Figs. 1, 3, and 5 (b) 
locate the known redshifts of 15 observed GRBs so far. The redshift values are quoted from 
http://cossc.gsfc.nasa.gov/batse/counterparts/GRB_table.html. We compare the observed 
redshift values only with the long subsample because all GRBs with the known redshifts 
belong to the long bursts. It is interesting to note that, if we assume a = 1.0, a large 
fraction of GRBs are distributed at high redshifts, provided that GRBs are assumed to 
follow the SFR, regardless of specific beaming models. However, in the case of a — 2.0, z max 
is significantly reduced, though it still remains quite large. The fraction for / >z / at z' = 3.42 
for a cylindrical-beam case is ~ 75 % for long and total samples when a — 1.0. This is 
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too large a value given the fact that the observed GRBs with such high redshifts account 
for < 10 % of the total bursts. This fraction become < 50 % when a = 2.0. It is also 
interesting to note that, though their (V/V max ) value is close to a Euclidean value, most of 
the short bursts are distributed at high redshifts in the case of a — 1.0. When we assume 
that the long bursts are uniformly distributed in space, the z max estimate for the broad 
beaming cannot explain the largest observed redshift in both cases of a = 1.0 and 2.0. 

Our study is based on an assumption that all GRBs have the same intrinsic luminosity 
and they are all beamed. Because the {V/V max } test is not sufficiently sensitive to the 
observed data, different luminosity functions are essentially indistinguishable. If there is 
an intrinsic luminosity function and the degree of beaming of GRBs is moderate, we may 
obtain results similar to what we have presented in this study. If all GRBs are highly 
beamed and their sources follow a SFR-like distribution, the maximum detectable redshift 
becomes very large, even in the case that the observed (V/V max ) is sufficiently close to 
the Euclidean value, i.e., 0.5. Therefore, it is difficult to rule out a possibility that the 
apparent Euclidean value of (V/V max ) may be due to the luminosity function, for instance, 
induced by beaming. In other words, the so-called Euclidean value may have nothing to do 
with the Euclidean distribution. If the total burst sample is a mixture of broad (or nearly 
non-beamed) beams and narrow beams, the exact fraction of the strongly beamed bursts 
may affect the burst rate estimate significantly. This is beyond the scope of this paper. 

For a given Lorentz factor, the range of the derived beamed luminosity distribution 
becomes broader as the beam opening angle decreases. We obtain log L max /L min ~ 13 in 
the cylindrical-beam case. However, in the case of broad beaming, log L max /L min becomes 
substantially smaller. This is only an effective range because the cut-off for the minimum 
luminosity given by the conic beam is rather loosely constrained. It is obvious that the 
luminosity range decreases when the opening angle increases for a fixed Lorentz factor. 
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This is simply explained as follows: Assuming that we observe a beamed emission with an 
opening angle A9, if we define ip as an angle between the line of sight and the symmetry 
axis of the conic beam, we can observe the beamed emission only within a certain range 
of ■?/>, Aip. For instance, when I/7 <S A9, Aip is approximately A9. For a given Lorentz 
factor, Aip increases as the opening angle of cone becomes larger. The apparent luminosity 
range from a broad beam would decrease gradually with ip, and in the extreme case, the 
isotropic radiation would have constant luminosity regardless of ip. However, when A9 is 
very small, the apparent luminosity rapidly decreases with ip, and the derived luminosity 
range is narrower than a broad-beam case. 



5. Conclusion 

We have demonstrated that the beaming- induced luminosity function may explain the 
statistics of the observed GRBs. In the case of the cylindrically beamed luminosity function, 
we determined the intrinsic luminosity L int by the (V/V max ) test for each subsample. In 
all cases, the intrinsic luminosity is much smaller than that of non-beamed luminosity 
functions, i.e. L ~ iq 52-53 erg/sec. The maximum luminosity obtained from our luminosity 
functions, in both the cylindrical-beam and the conic-beam cases, is L mSLX ~ iq 50 " 51 erg/sec. 
This is compatible with the required value for isotropic radiation. The ratio between the 
maximum luminosity for the long bursts and that for the short bursts -^ m ax,ion g /-^max,short 
is ~ 7 for the narrow beam and ~ 3 for the broad beam, respectively. This implies that 
these two subgroups may have two different intrinsic luminosities, and hence probably two 
different origin. The obtained z max in the cylindrical-beam case are ~ 10, 14,3 (a = 1.0) 
and ~ 4.5,5.5, 1.6 (a = 2.0) for the total selected sample, and long and short subgroup, 
respectively. The obtained z mSLX in the broader beam case (i.e. A9 = 3°.0) are ~ 4,6,2 
(a = 1.0) and ~ 2.2,2.6, 1.2 (a = 2.0) for the total selected sample, and long and short 
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subgroup, respectively. The maximum known redshift we consider, z = 3.42, is estimated 
for GRB971214. However, even a larger value of z < 3.9 is suggested for GRB980329 
(http://cossc.gsfc.nasa.gov/batse/counterparts/GRB_table.html). Without any detailed 
information on the redshift distribution of GRBs, the maximum redshift gives a rather 
strong constraint on the ratio between the Lorentz factor and the opening angle A9. For 
the conic beam case, the luminosity function derived for the narrow opening angle seems to 
fit the (V/V max ) curve better. 

From this simple analysis we have carried out, we have shown that the beaming-induced 
luminosity function may account for the basic statistical properties of the observed GRBs. 
The apparent Euclidean distribution of GRBs may be an indication of the presence of the 
GRB luminosity function. 

Although it is the beyond the scope of this work, it is interesting to point out that there 
could be a potentially observable correlation between spectral hardness and luminosity 
since the most Doppler-boosted burst would probe the hardest part of the spectrum. Any 
wide scatters in intrinsic spectral shapes would make the detection of this correlation 
complicated. 

We thank H. Kim and K. Kwak for useful discussions. IY is supported in part by the 
KRF grant No. 1998-001-D00365. We appreciate the referee, Ralph A. M. J. Wijers, for his 
suggestions and some useful information. 
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Fig. 1. — Two statistical properties of three sample sets for the cylindrical-beam case are 
shown. Data sets are indicated at the upper right corners in each panel. In the left panels, we 
plot the (V/Vmax) curve as a function of number of bursts. Dotted lines represent the observed 
(V/V max ) with ±3er bounds. Solid lines and dashed lines show theoretical results, the former 
result from the SFR-motivated number density distribution for the spatial distribution of 
GRBs and the latter the uniform distribution. In the right panels, we show the f >z > in terms 
of z' for each sample set. The y-axis is the fraction of GRBs which would be located farther 
than a certain redshift z'. Solid lines and dashed lines are same as in the left panels. Filled 
squares indicate the estimated redshifts of 15 observed GRBs. In all panels, the photon index 
a is assumed to be 1.0 (thin lines) and 2.0 (thick lines). Note that the z max is significantly 
reduced when a = 2.0. 

Fig. 2. — Plots simlar to Fig. 1. Results for the total sample for different opening angles 
in the conic beam model. Opening angles are indicated at the upper right corners in each 
panel. In all panels, the photon index a is assumed to be 1.0 (thin lines) and 2.0 (thick 
lines). 

Fig. 3. — Same plots as Fig. 2. for the long bursts. In all panels, the photon index a is 
assumed to be 1.0 (thin lines) and 2.0 (thick lines). 

Fig. 4. — Same plots as Fig. 2 for the short bursts. In all panels, the photon index a is 
assumed to be 1.0 (thin lines) and 2.0 (thick lines). 

Fig. 5. — Same plots as Fig. 2 for the conic-beam case for the total sample. We adopt a 
bi-Gaussian function as a number density function of GRB sources, n(z), which is different 
from that used in Figs. 1-4. It gradually decreases at high redshifts. The photon index a is 
assumed to be 1.0. 
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